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Introduction 


GMAO is close to having a working 4DVAR system and has 
developed a linearized version of GEOS-5. 

This talk outlines a series of improvements made to the linearized 
dynamics, physics and trajectory. 

Of particular interest is the development of linearized cloud 
microphysics, which provides the framework for ‘all-sky’ data 
assimilation. 
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• Introduction 

• Updated Dynamics 

• Updated Boundary Laye 

• Advection 

• Moist Physics 

• Radiation 

• Gaussian Quadrature 

• Sensitivity Studies 


Outline 
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Metrics for improvement 

The linear model is validated by comparing correlations, RMSE 
and RMS between the nonlinear perturbation trajectory, 

m (x + <5x) — m (x) , 

and the tangent linear model perturbation trajectory, 

M5x, 



Time (UTC) 
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Updated Dynamics (Jong Kim) 


Noticed that with increased resolution the linear model has issues 
at the jet levels. 

• Update to sponge levels to make more like NL model. 

• Reduced ‘KORD’ in remapping (more stable and linear). 

• Added option to apply Shapiro filter at the jet levels 
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Updated Dynamics - RMSD 






X ltf 6 


6/37 


Height (hPa) 


Updated Dynamics - RMS 
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Height (hPa) 


Updated Boundary Layer 


Originally a simple Louis type BL scheme was used in the linear 
model to stabilize the dynamics. 

• Large correlation reduction near the surface for winds and 
temperature. 

• Potential to upset the cloud scheme. 

(r) 

Updated to use the full K y m h that is used in the nonlinear model. 

• Perturbations K' m h still neglected to avoid steep gradients near 
stable-unstable (/?,- = 0) transition. 
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Updated Boundary Layer - Correlations 


(a) u' (ms ! ) 



(b) v' (ms ') 



Correlation 


(c) t; (K) 



Correlation 


(d) q' (kgkg l ) 
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Linearized Advection 


Another aspect of the model that interacts strongly with 
prognostic clouds is the advection. 

Tracers in linearized GEOS-5 advected with schemes used 
nonlinear model. 

• PPM of Colei la and Woodward (1984)/Van Leer (1977) 

• Uses Huynh 1996 second constraint and Lin 2004 RMC. 

These schemes have the potential for nonlinearity. 


Passive cloud test 


To test the tangent linear behavior of the advection we set up a 
passive cloud experiment. 

• Only perturbed variables are clouds. 

• Perturbation is infinitesimal and has structure of field itself. 

• Cloud uncoupled from other fields (GOCART and Radiation off). 

• Wind is constant in time ( u ' — 0) so continuous advection is 
linear. 
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Passive cloud 


correlations 
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NLM ql (PPM Lin Scheme) 
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Passive cloud 


correlations 
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NLM ql (PPM Lin Scheme) 
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Advection ID Case Study 

We don't really want to put up with either of these situations so 
we’re investigating further... 

[ K = u (r)W 

dt dx ’ 

u = 1, x = (0, 1], N = 64, Ax = 1/A/, At = 0.1 / N. Three test 
cases. 


0.9 
0.8 
0.7 
0.6 
cr 0.5 
0.4 
0.3 
0.2 
0.1 
0 

0 0.2 0.4 0.6 0.8 1 

X 

Perturbations are initial profiles multiplied by 1 x 10 -4 


Initial profiles 



16/37 


Advection ID Case Study 


3rd order upwind 


PPM With Lin 
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In Matrix form, 


Advection ID Case Study 




5 


Seek solution q ~ exp( — /At) 


-Aq' = Aq'. 


Test case 2 PPM scheme: 


Eigenvalues 
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Eigenvector 1 
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Linearized cloud scheme 


Linear cloud scheme to accompany recently implemented linear 
RAS scheme Holdaway et al (2014). 

• Based on Bacmeister scheme used in nonlinear model. 

• Single moment microphysics (mass mixing ratio). 

• Prognostic treatment of clouds. 

• Linear model variables q^, q' : and C^ N . 

• Cloud scheme variables q\ LS , q\ AN , q\ LS , q[ AN , C' LS and C' AN 

• Trajectory q [ r $ , q { AN and C^. 

• Sources using RAS &i PDF approach, evaporation, sublimation 
autoconversion and falling and reevaporation of precipitation. 


Problems in large scale cloud source 


The cloud scheme uses many nonlinear functions with 
discontinuities and steep gradients. 

Very large perturbations occur around the calculation of cloud 
fraction when using PDF approach. We develop two approaches to 
dealing with this: 

• Perturbation model approach for use in 4DVAR 

• Filtering + perturbation model for use in observation impacts 
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Large scale cloud source perturbation model 
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Latitude (°) Latitude (°) Latitude (°) 


6-hour cloud perturbations 


(a) q ( Nonlinear Perturbation at 850hPa 
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(a) q. Nonlinear Perturbation at 250hPa 
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(b) qj Tangent Linear Perturbation (Moist) 
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6-hour correlations and RMS 


Correlations 


RMS 


(a) q 1 (b) q. 




(a) qj (b) q. 
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Height (hPa) 


Filtering 

The perturbation model does not perform well for longer than a 
12-hour integration. Large perturbation growth vs. way too 
diffusive. 


Filtering developed based on examining the Jacobian of just the 
large scale source routine. 
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Filter based on magnitudes of the elements and eigenvalues. 
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xe/Je xe/oe ' ie/re 


Jacobian correlations 
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24-hour correlations - Global 
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Height (hPa) 


24-hour correlations - Tropics 
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Height (hPa) 


Radiation 


Solar and infrared radiation is now an option in the linear model. 

Derived exactly from nonlinear code so slow, for the adjoint it’s 
more expensive than the dynamics! 

• Connected to clouds mass and fraction but not precipitation. 

• Ozone is only active chemistry, CO 2 global constant. 

• Option to include dust but 5 bins so really slow all round. 

...nothing to see here... 
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Gaussian quadrature - Rahul Mahajan 


Observation impacts at GMAO are computed using two adjoint 
sensitivities. Linearized around forecasts initialized with the 
background and analysis. 

Using Guassian quadrature the adjoint is linearized around an 
average of the two. 

• Improves correlations for 24-hour forecasts. 

• Extends linear validity to 48-hours so we can do longer 
observation impacts. 
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24-hour correlations w/ Gaussian quadrature 
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Height (hPa) 


Adjoint sensitivity case studies 


Interested in using the adjoint to analyse: 

• Predictability of storm track and intensification (Brett Hoover) 

• The role of dust in TC development (NAMMA period) 

• Sudden stratospheric warming (Jan. 2013 and Jan. 2009) 
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Dust Case Study 


The development of linearized radiation enables the investigation 
of sensitivity to aerosols. 

• NAMMA period mid August to mid September 2006 

• Significant dust outbreaks combined with TC cyclogenesis 

• Significant debate over the role that dust plays. 

• Particular interest in hurricane Helene 

Start by examining the following: 

• August 27th OOz to August 28th OOz. 

• Forecasts and adjoint at C360 (25km) with linear radiation. 

• Currently no linear GOCART. 
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Sea level pressure (Pa) 
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Conclusions 


• Improvements to the current dynamics and boundary layer 
scheme. 

• Significant speed up of convection by moving from Taf to 
Tapenade. 

• Linear clouds in and working. Next step is all-sky assimilation 
(Min-Jeong Kim @ AGU). 

• Radiation with dust. 

• Have the confidence to look at some sensitivity studies. 

• ...But question remains over advection, would like to improve. 
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Questions 
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